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Abstract We show that an integro-differential equation 
model for pulse propagation in optical transmission lines with 
dispersion management, is integrable at the leading nonlinear 
order. This equation can be transformed into the nonlinear 
Schroedinger equation by a near-identity canonical transfor- 
mation for the case of weak dispersion. We also derive the 
next order (nonintegrable) correction. 



I. INTRODUCTION AND MAIN IDEAS 

The impressive progress in high-bit-rate optical fiber 
communications achieved during the last few years was 
primarily due to the invention of the dispersion man- 
agement (DM) technique Q, which has become the key 
technological component of fiber telecommunications. A 
recent example of just how successful this technique is 
the 1.02 Tbit/s data transmission achieved in standard 
optical fiber with dispersion management over 1000km 
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The main factor that limits the bit-rate is pulse- 
broadening due to the chromatic dispersion of optical 
fiber. This broadening is characterized by the dispersion 
length, Zdis ~ (d x (BR) 2 ) -1 . Here, d is the fiber chro- 
matic dispersion and BR is the bit-rate. The dispersion 
length Zdis is the distance at which the pulse-width ap- 
proximately doubles due to dispersive broadening. This 
distance decreases as a quadratic function of the bit-rate. 

The basic idea of the DM technique is to compensate 
the fiber chromatic dispersion by periodically incorporat- 
ing an additional element, such as a dispersion compen- 
sating fiber or fiber chirped gratings, whose sign of the 
chromatic dispersion coefficient is opposite to that of the 
main optical fiber. (See Fig.). 



► Optical Amplifier 

~ Compensating Fiber 

-2*2 Transmission Fiber 



FIG. 1. A schematic of the dispersion management of an 
optical fiber link is presented on the upper part of the figure. 
The dependence of fiber chromatic dispersion Dasa function 
of the distance z is presented below. 

Modern optical transmission systems must satisfy very 
strict requirements for bit-error-rate (BER ~ 10~ 12 to 
10 -15 ). Therefore, the pulse amplitude should be large 
enough so that it can be effectively detectable. As a con- 
sequence, the Kerr nonlincarity of the fiber refractive in- 
dex n = rio + al (uq - linear part of the refractive index, 
/ - pulse intensity, and a - coefficient of Kerr nonlin- 
earity) should be taken into account. The spectrum of 
an optical pulse with characteristic power Pq will expe- 
rience noticeable nonlinear distortion at distances grater 
than the characteristic nonlinear length Z n \ = (cvPo) -1 - 
A natural idea is to control the nonlinear spectral dis- 
tortion by properly choosing the value of the residual 
dispersion over the period of the fiber link. 

The dispersion management technique is currently the 
focus of intensive experimental, numerical and theoretical 
research p|-pd| because of its excellent practical perfor- 
mance. Optical pulses in DM fiber links exhibit unex- 
pected soliton-like properties, such as high stability and 
elastic interaction, despite the fact that the governing 
equation does not belong any known class of integrable 
equations. 

In this paper, we demonstrate that the leading-order 
equation describing the slow dynamics [Q of optical 
pulses in such systems is close to integrable for the case of 
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weak dispersion management. We also present the next 
order (nonintegrable) correction. This correction is small 
for weak dispersion management. 

The problem of weak dispersion management was con- 
sidered in |Q for the first time as a theoretical exam- 
ple of periodic variation of the dispersion coefficient in 
the model of transmission link. The power of Lie trans- 
form technique proposed by authors earlier in Ju^,[w| was 
also shown in |L2j . This paper was published long before 
dispersion compensation became popular and before its 
strength was proven experimentally. The main result of 
this paper was that pulse propagation in fiber links with 
weak variation of dispersion can be described in leading 
order by an unperturbed nonlinear Schroedinger equa- 
tion. Using a formalism which is common in weak turbu- 
lence theory [^5), we derived similar results and investi- 
gated its validity by calculating the higher order correc- 
tions. 



II. BASIC EQUATIONS 

The relevant equation for the electric field envelope 
E(z,t) is the Nonlinear Schroedinger equation (NLS) 
with periodically varying dispersion and an external force 
representing fiber losses and amplification: 

iE z + \^d{z)E tt + \E\ 2 E = iR'{z)E (1) 



2 Zdi 
R'{z)E = Z nl 



~7 + r^S(z - z k ) 



k=l 



(2) 



Here, 7 describes the fiber losses. The amplifiers that 
compensate the fiber losses are placed periodically at 
points Zfe and separated by distances Z a ~ f /j. The coef- 
ficient of amplification r is chosen to exactly compensate 
energy losses after passing the amplification distance Z a . 
The function d(z) is defined by 



d{z) 



f, in the transmission fiber, 



compl dtransi 



in the compensating fiber 



Without loss of generality, we consider the case when 
the period of amplification is equal to the period of com- 
pensation. The theory can be trivially extended to the 
general case. 

In practice, Z n \/ Z^ s ^> I and jZ n \ 3> I. Therefore, 
on the scale of the amplification distance Z a , the pulse 
dynamics is practically linear. On the other hand, the 
linear evolution of the signal preserves the Fourier spec- 
trum of the pulse. One can consider optical pulse prop- 
agation over one period of a fiber link to be a mapping 



of the input pulse into the output pulse. To preserve the 
bit pattern, this mapping must be both stationary and 
stable. Stationarity can be achieved by a proper choice 
of the amplification coefficient of the amplification and 
by complete compensation of dispersion. In this case, 
stationarity will be achieved for any pulse shape. Stabil- 
ity, on the other hand, is more complicated, and can be 
investigated through the slow nonlinear pulse dynamics. 

In order to consider slow nonlinear pulse dynamics it 
is therefore natural to transform the equations into the 
Fourier representation, and study the slow dependence 
of the spectrum on the coordinate z. One can see the 
separation of scales in the equation (|f]) : E z is dominated 
by the RHS and dispersion term. Let us remove the 
external force (RHS) from |l]) by the transformation E = 
q(z, t) exp (i?(z)), and we obtain the equation 



i( lz + \-^-d{z)q tt + c(z)\q\ 2 q = 0, 

where c{z) — exp (2R(z)). 

Let us decompose d(z) into the sum 

d{z) = (d(z)) + d(z), (d(z)) = 0, 



(3) 



(4) 



where 



</(*)> = y I f ^ dz - 



Following §], we represent q(t, z) as 



q(z,t) = / duja UJ (z) exp (iuit + ioj 



2 1 Znl 



2 Zdis Jo 



d(z)dz), 



and compute the equation for the slow dynamics of the 
pulse's Fourier spectrum in Hamiltonian form 



l-+«77 = 0. 

oz da* 



(5) 



with the Hamiltonian 

H = [ 7:77— < d(z) > ^la^duj + [ dujiduj2duj3duj4 
J 2 Zdis J 

xF ^4 S ( LJ i ^4)a* Ul a* 2 a U3 a Ui , (6) 

where 

F^! = c(z) exp + u 2 - - f d(f)de). 

^ ^dis Jo 
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Note that the kernel depends only on the com- 

posed variable g = (wf +w| — w| —uj\), which, as we will 
see later, leads to important consequences. 

The leading-order approximation for the slow nonlin- 
ear dynamics can be obtained by averaging (^|) over Z a . 
The resulting equation is valid if the pulse spectrum does 
not change appreciably over the amplification distance. 
The next order correction to the averaged equation was 
calculated in jl6| by using a Lie transform jl3|]l4| ]. An 
alternative method for calculating these corrections was 
proposed in jlTj . 

In the simple case of communication through one fre- 
quency channel, with the additional assumption that 
pulses do not interact with each other, one can use the 
quasi self-similar character of the single-pulse behavior 
p9[ , which was proposed for the first time by Talanov 
(lens transformation) ptJj. As a result, the slow dy- 
namics of a single pulse is described by the Nonlinear 
Schroedinger equation with a quadratic potential. This 
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approach was extensively discussed in Q , |l9| 
allows one to explicitly take into account the evo 
of the single pulse spectra within the system period. An- 
other advantage of the quasi self-similarity is that it also 
holds for very small values of the residual dispersion. It 
was shown in |2^] that stable single pulse propaga- 
tion is possible for zero or even negative residual disper- 
sion for special configurations of the dispersion map. 

However Talanov's method cannot be applied to the 
case of several interacting pulses within one frequency 
channel, nor can it describe the evolution of an optical 
pulse propagating in a system that utilizes the wave- 
length division multiplexing (WDM) technique. The 
WDM technique exploits several frequency channels for 
data transmission. Optical pulse streams in systems that 
use this technique propagate in every channel with their 
own velocities and interact through the common refrac- 
tive index. 

In the fiber links with dispersion management, the 
nonlinear effects are smaller than the dispersive effects, 
therefore the Fourier spectrum of a pulse does not change 
much between amplifiers. Hence, equation (^) can be re- 
placed by the equation averaged between the amplifiers. 
In the corresponding averaged Hamiltonian, F£^ [ s re _ 

>, so that this Hamiltonian 



placed by 
becomes 



-<- F<*>lV2 



TL = I -77— < d(z) > cj 2 |a w | 2 (iti) + j dioxdw^dLOzdiOA 



here 



rpuJ±UJ2 — 



-a Jo 



dzc(z) 



dis Jo 



This averaging is equivalent to assuming that higher or- 
der corrections in the Lie expansion are small compared 
to the leading order. 



III. INTEGRABILITY AT THE LEADING 
ORDER 

In our further analysis we will utilize ideas that are 



well developed in the theory of weak turbulence 15 1. We 
rewrite the Hamiltonian (^) as 



n 



QuCL^a* du 



a* t a Z 2 % Sui+ui-us-ui duJidu) 2 du 3 du)4 , 



(8) 



where we define the leading order quadratic dispersion 
relation as 

Q u =f3<d{z)>u\ (3=1^, (9) 

^ ^dis 

Then a w satisfies the canonical equation of motion: 

If T^^2 is a constant independent of the w's, then equa- 
tions (PHlO|) correspond to the usual NLS equation. With 
T^ 4 2 given by (0) || , these equations describe the aver- 
aged, slow-time evolution of an optical pulse propagating 
in a fiber-optical link with periodic amplifiers, damping, 
and dispersion compensators. 

To find the leading-order approximation to (|^-^0|), we 
follow [fL5 23 25|] and perform the transformation 



&uj + J B LU2UJ3 b LUl buj 2 b L j 3 S UJ -^ U j 1 — cj 2 — uj 3 + 

higher order interactions. (11) 



The transformation ( pL l| ) is canonical up to and includ- 
ing the cubic terms if the coefficients B^d, satisfy the 
symmetry condition fig] 



Substitution of ( |11[ ) into the Hamiltonian (|8|) transforms 
(I) into 
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H 



Qu>b UJ b*dw- 



K t K 2 bu 3 bu>i 5 m +UJ2 - ul3 - UJl duidu 2 du; 3 dui4 
+ higher order interactions. (12) 

The main mechanism of wave interaction is four-wave 
scattering, which occurs when ui's and Q^'s satisfy the 
following resonant conditions 



Qu "i" Qui — 



(13) 



The set of oj's that satisfies these conditions is called the 
resonant manifold. The form of the transformed Hamil- 
tonian (|l2|) explains the importance of the notion of a res- 
onant manifold. Indeed, at the points where (|l^) holds, 
in general, one cannot choose the arbitrary coefficients 
B %%t>\ to exclude the terms a* 1 a* 2 a^ 3 a„ 4 from the 

Hamiltonian because of "small resonant denominators" . 

For the case of quadratic dispersion (|j|) , equations ( |T3] ) 
have only trivial solutions 



LO2 

L0 3 



UJ 3 
UJ2 



(J, 
U). 



(14) 



T(9)=T(9 = 0)+ 9 -+Y: 

y k=2 



k fjrpuiUli 

~kl dg 



(17) 



T(g = 0) = To, i.e. a constant independent 
of uji . The transformed Hamiltonian acquires the form 



Here T" 1 " 1 



H 



T b*b* b^b^A 



Ul U U>1 U ^2 "U+Ul — UJ2 —U! 3 



Qubub*dio 



higher order interactions, 



(18) 



which is the Hamiltonian of the NLS equation to lead- 
ing order. Thus equation ( p^| ) with Hamiltonian (^|) and 
quadratic dispersion relation (^|) is integrable at the lead- 
ing nonlinear order. 

The condition of quasi-identity of transformation jll| ) 
reads as 



Jj" 1 " 2 < 1 



(19) 



Let us estimate the in (|1 



For the case of a loss- 
less piece-wise constant dispersion map with alternating 
pieces of fibers of length li with dispersion < d > +di and 
of length I2 with dispersion < d > +c?2, (Wi + hd 2 = 0), 
T" 1 " 2 acquires the form 



In this case the value of T^) on the resonant mani- 
fold (|13|),(|14|) is equal to a constant independent of u/s: 
= T . We will show below that the leading nonlin- 
ear approximation the equation ( |To| ) with Hamiltonian 
(H) and quadratic dispersion relation (^|) is the NLS, and 
hence it is integrable. Thus, we can integrate ([!(]) at 
the leading nonlinear order. 

For the case of interest, the interaction matrix element 
is a function of g = (u>\ + uj\ — uj\ — ): 



T(g). 



Now we choose the canonical transformation kernel as 



or equivalently, 



CJ3LJ4 UJ\U)\ 



T(g) - T(0) 



(15) 



(16) 



Note that this expression does not have a "small denom- 
inator" problem. It is not singular on trivial resonances, 
i.e. it is not singular for g = 0. This fact is obvious from 
an expansion of m powers of g: 



rpLJ\U2 



sin(tfp) 



(20) 



where S is the strength of a dispersion map. Thus condi- 
tion ( |Tg| ) becomes 



S 2 « (d), 



(21) 



so that the transformation ( |ll| ) is quasi-identical for the 
so called weak dispersion maps, i.e. for the cases of small 
variation of the dispersion on top of big dispersion. 

In this paper, we only study the sixth-order interaction 
in the averaged equations. V.E. Zakharov |l7j recently 
pointed out that a canonical transformation similar to 
( |TT| ) can also be used to bring the original nonaveraged 
problem (^|) into a form analogous to that developed 
here for the averaged equation, namely, the Nonlinear 
Schroedinger Equation with higher-order corrections. We 
believe that the present work is an important stepping 
stone in the proper development of a more sophisticated 
transformation, which will give us sufficient understand- 
ing to attack this much harder problem. 
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IV. SIX- WAVE INTERACTIONS ON A 
RESONANT SURFACE 



In this section, we calculate the sixth-order interaction 
matrix element for the Hamiltonian (|l8|). The Hamil- 
tonian (|l8| ) is the averaged Hamiltonian, and, including 
sixth-order terms, has the form: 



H = / dtoQJb^ 



T Q dwdu: 1 dhJ' i dhJ i bl ]1 bl ]2 b UJ3 b uli 



dwduj i dtv 2 doj 3 doj 4 du 5 T^^^ 



h*h* h* h h h A w +"i+ w 2 



(22) 



There are two ways of calculating the sixth-order in- 
teraction matrix element T^.f, };f 2 for the Hamiltonian 
(|22|) . The first way is to calculate the higher order terms 
in the canonical transformation (|TT|) . The second way, 
which is simpler and clearer, is to compute r,""?"? g in- 
volving standard Wyld diagram formalism |18j which is 
a generalization of the Feynman diagram technique for 
classical noneq uilibrium systems. 

Following |2(| the matrix element T££*%* can be rep- 
resented as sum of all topologically different irreduceable 
diagrams with no internal loops. Below we present, an 
example of such a diagram. The structure of this diagram 
can be understood as follows: This diagram consists of 
two vortices A and B. Each of these vortices forms a di- 
agram which corresponds to four wave mixing. For com- 
puting of six order interaction term in the perturbation 
theory one should take two four-wave vortesces, connect 
them using all topologically possible configurations and 
then summ all of them. The number of such nonequva- 
lent configurations is nine. The diagram presented below 
is one one of these nine terms. 




The analytical expression corresponding to this dia- 
gram is 



^OJl, 0^2+^3— (jJ6 r J lUJ 2 iJ 3 



Qu)2 ~t~ QuJ3 QuJq Qw2+^3-W6 

After summing up all 9 diagrams, the sixth order ma- 

i=9 

/ In. ! - 

(23) 



trix element T^J.fJ can be represented as 



rriUiU2U3 



£ 

i=l 



\T(m) 



where the nine variables gi 1 
following expressions 



1, .., 9 are given by the 



{-2(w 4 


- wi)(w 4 


- w 2 ), -2(cj 4 


- wi)(w 4 


-w 3 ), 


-2(cj 4 


- UJ 2 ){uJi 


- w 3 ), -2(w 5 


- wi)(w 5 


-w 2 ), 


-2(cj 5 


- <*>i)(w 5 


- uj 3 ), -2(uj 5 


- u 2 )(u> 5 


-w 3 ), 


-2(w 6 


- wi)(w 6 


- uj 2 ), -2(uj 6 


- Wi)(w6 


- w 3 ), 






-2(w 6 - 


- W 2 )(u) 6 - 


- w 3 )}. 



The resulting expression can be viewed as the measure 
of the "distance" of the perturbed system from its inte- 
grable analog. Indeed direct calculations show that for 
the case of NLS, when T^ 1 '^ 2 = const, this sixth order 
matrix element is identically equal to zero on the reso- 
nant manifold 



Qu ~l~ Qui ~l~ Quj^ — Qujs ~t~ Qui Qujq 

UJ + Ull + Ul 2 = UJ3 + OJ4 + 0J5- 



In order to bring equation 



(24) 

J) as close to the integrable 



limit as possible, we must minimize T££J£jj} over all the 
tunable system parameters. 

For the case of a lossless piece- wise constant dispersion 
map, when the fourth order matrix element is given by 
(pp|), we find (using Mathematica 3.0), that the six wave 
interaction matrix element 

T^ 5 ^S 2 /(d), 

on the 6- wave resonant manifold (pi]). Therefore the 
above perturbation technique is valid for the case of a 
weak dispersion map. In leading order, the slow pulse 
evolution in optical fiber is governed by the integrable 
NLS equation. 
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V. CONCLUSION 

We showed that the leading-order equation describ- 
ing pulse dynamics in optical fiber links with a weak 
dispersion map is close to integrable. Using canonical 
transformations, we found a nonintegrablc higher order 
correction to the integrable part of the model equation 
describing the slow evolution of optical pulses. 
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